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Abstract 

An extension to a high-order model for the direct isotropic/smectic- 
A liquid crystalline phase transition was derived to take into account 
thermal effects including anisotropic thermal diffusion and latent heat 
of phase-ordering. Multi-scale multi-transport simulations of the non- 
isothermal model were compared to isothermal simulation, showing that 
the presented model extension corrects the standard Landau-de Gennes 
prediction from constant growth to diffusion-limited growth, under shal- 
low quench/undercooling conditions. Non-isothermal simulations, where 
meta-stable nematic pre-ordering precedes smectic-A growth, were also 
conducted and novel non-monotonic phase-transformation kinetics ob- 
served. 

1 INTRODUCTION 

The kinetics of phase transformations [f ] is a fundamental subject in mate- 
rial and interfacial science that has widespread impact on material man- 
ufacturing and use. Three fundamental kinetic phenomena associated 
with phase transformations are microstructure, growth rate, and shape 
evolution [2] (see [I]). The growth processes associated with phase trans- 
formations vary, depending on driving forces, such as bulk free energy 
minimization and reduction of interfacial area. Growth laws in diffusive 
and non-diffusive transformations typically differ in that the former has 
a conserved order parameter [1,2]. The differing kinetics of non-diffusive 
and diffusive phase transitions are manifested in their growth laws I oc t n , 
where I is the characteristic length of the growing domain. Non-diffusive 
phase transitions exhibit constant growth where I scales linearly (n = 1). 
Diffusive phase transitions exhibit diffusion-limited growth, where n = | . 
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Figure 1 : Schematic showing the inter-relationships between the phase-ordering 
processes. Adapted from Figure 1 of rcf. [35] 

In the case of liquid crystals, both (mass and thermal) diffusive and 
(phase-ordering) non-diffusive transformation dynamics are present simul- 
taneously. The study of diffusive transformations in liquid crystalline 
materials has mainly focused on mass diffusion. Two main areas of this 
research involve liquid crystalline materials where either a small amount of 
(undesired) impurity is present or where a composite material is desired, 
such as polymer dispersed liquid crystals [3] (PDLCs). In the former case, 
much work has focused on the study of diffusion-driven growth instabil- 
ities [4,5] (see refs. [6,7] and Chapter B.IX of ref. [8]) in nematic liquid 
crystals both through experimental and theoretical approaches [9-16]. 
Higher order mesophases such as smectic and columnar liquid crystals 
have been studied, in this context, to a lesser extent [17,18]. Much of 
this interest stems from the fact that liquid crystal transformations occur 
on experimentally accessible timescales and are inherently anisotropic. In 
the latter case, the study of PDLC materials [3] has focused on differ- 
ent mechanisms of phase-separation [19-28] (polymerization, solvent, and 
thermally induced), where growth morphology of liquid crystal-rich do- 
mains and their texture are studied in order to control functional physical 
properties, such as electro-optical response. Subsequently, PDLCs com- 
posed of higher order liquid crystals have also been studied to a lesser 
extent [29-31]. 

On the other hand, thermal diffusion in liquid crystalline phase tran- 
sitions has typically been neglected due to the relatively low latent heat 
contributions of (nematic) liquid crystalline transformations and the per- 
vasive use of thin film geometries [7] (where evolved latent heat escapes 
in the vertical dimension). Focusing on pure/single-phase liquid crystals, 
where impurity concentrations are below the saturation limit, evidence has 
shown that latent heat effects can also result in diffusive dynamics [32]. 
Thus heat diffusion must be taken into account to reproduce experimen- 
tally observed growth laws under these conditions [33,34]. As the study 
of these non-isothermal effects progresses [32,35,36], it is becoming clear 
that they play a non-negligible role in the complex kinetics and dynamics 
of liquid crystalline phase transitions. 

As previously mentioned, a relatively large amount of effort in the 
study of liquid crystals has focused on those which exhibit only some de- 
gree of orientational order, or nematics. Higher-order mesophases that 
exhibit some degree of positional order, in addition to orientational, have 
been less studied. The most simple of these higher order liquid crystals 
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is the smectic-A mesophase, which exhibits lamellar translational order- 
ing, in addition to the orientational ordering of nematics. Recently, an 
increasing amount of interest in this mesophase, in particular of materials 
exhibiting a direct isotropic/smectic-A (disordered/ordered) transition, 
has resulted in many experimental and theoretical results. Nonetheless, 
the understanding of this mesophase is in a nascent stage. Much of this 
is due to the time and length scales at which the structures and dynamics 
occur being on the nano-scale. 

The additional positional order of smectics, where the underlying pri- 
mary order is orientational, allows the possibility of dual-order kinetics 
and dynamics. This has been experimentally observed in a smectic poly- 
mer liquid crystal which exhibits extremely slow dynamics [37]. The 
presence of a transient meta-stable nematic phase preceding the stable 
smectic phase in this polymer liquid crystal was observed under certain 
quench conditions [38]. The general phenomena of meta-stable states 
in phase transitions involving dual non-conserved order was theoretically 
shown over a decade ago [39], but this generalized approach is not suit- 
able for liquid crystal phase transitions |T]). A first-approximation of the 
experimental system [37], taking into account shape kinetics/growth ki- 
netics/texturing (see [I]), has been modeled [40] shedding light on some of 
the nano-scale phenomena resulting in these experimental observations. 

In this context, the kinetics of liquid crystal phase transitions, in par- 
ticular smectics, becomes more complex and diverse than currently re- 
garded where coupled non-diffusive (phase-ordering), diffusive (mass/thermal), 
and pre-ordering transformation dynamics are present. The general ob- 
jective of this work is to contribute to the fundamental understanding 
of mesophase formation under simultaneous orientational and positional 
symmetry-breaking and under non-negligible latent heat evolution. The 
specific objectives of this work are (1) to extend an existing high-order 
isothermal model for the direct isotropic/smectic-A transition to account 
for thermal effects and (2) to conduct a basic study of the phase transition 
landscape: 

1. Derive an energy balance to extend an existing high-order model 
for the isotropic/smectic-A phase transition [41,42] to account for 
non-isothermal effects including latent heat of phase-ordering and 
thermal diffusion. 

2. Simulate a one-dimensional growing smectic-A front, using the ex- 
tended model, to show that the non-isothermal extension corrects 
the standard Landau-de Gennes prediction of volume-driven growth 
kinetics (I oc t) to diffusion-limited growth kinetics (I oc t 1 ' 2 ) under 
shallow quench/undercooling conditions, consistent with experimen- 
tal observations [18,33,34]. 

3. Simulate a one-dimensional growing smectic-A front using the non- 
isothermal model under conditions where meta-stable nematic pre- 
ordering is observed [43]. 

This work neglects nucleation mechanisms [44], fluctuations [45,46], impu- 
rities [3,8], and convective flow [46] while taking into account energetically 
the inter-coupling between orientational/translational order and variation 
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Figure 2: schematics of the (a) isotropic, (b) nematic, and (c) smectic-A phases. 

of smectic layer spacing. This work is organized into effectively three sec- 
tions: background/model, energy balance derivation, and simulation. The 
background/model section provides a brief introduction to relevant types 
of liquid crystal phase-ordering and the Landau-de Gennes type model 
of Mukherjee, Pleiner, and Brand [41,42] used in this work. The deriva- 
tion of the thermal energy balance extension to this model is then derived. 
Finally, simulation results are presented and discussed, followed by a sum- 
mary of the conclusions. 



As mentioned above, liquid crystalline phases or mesophases are materi- 
als which exhibit partial orientational and/or translational order. They 
are composed of anisotropic molecules which can be disc-like (discotic) 
or rod-like (calamitic) in shape. Thermotropic liquid crystals are typi- 
cally pure-component compounds that exhibit mesophase ordering most 
greatly in response to temperature changes. Lyotropic liquid crystals are 
mixtures of mesogens (molecules which exhibit some form of liquid crys- 
tallinity), possibly with a solvent, that most greatly exhibit mesophase 
behavior in response to concentration changes. Effects of pressure and 
external fields also influence mesophase behavior. This work focuses on 
the study of calamitic thermotropic liquid crystals which exhibit a first- 
order mesophase transition. 

An unordered liquid, where there is neither orientational nor transla- 
tional order (apart from an average intermolecular separation distance) 
of the molecules, is referred to as isotropic. Liquid crystalline order 
involves partial orientational order (nematics) and, additionally, partial 
translational order (smectics and columnar mesophases). The simplest of 
the smectics is the smectic-A mesophase, which exhibits one-dimensional 
translational order in the direction of the preferred molecular orientational 
axis. It can be thought of as layers of two-dimensional fluids stacked upon 
each other. Schematic representation of these different types of ordering 
is shown in [2] 

Theoretical characterization of orientational and translational order 
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BACKGROUND AND MODEL 



2.1 



LIQUID CRYSTAL ORDER 
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of the smectic-A mesophase (see [2| is accomplished using order param- 
eters that adequately capture the physics involved. Partial orientational 
order of the nematic phase is characterized using a symmetric traceless 
quadrupolar tensor [41]: 



where n/m/1 are the eigenvectors of Q, which characterize the average 
molecular orientational axes, and S/P are scalars which characterize the 
extent to which the molecules conform to the average orientational axes 
[47-49]. Uniaxial order is characterized by S and n, which correspond 
to the maximum eigenvalue (and its corresponding eigenvector) of Q, 
S = f/^n- Biaxial order is characterized by P and m/l, which correspond 
to the lesser eigenvalues and eigenvectors, P = | (fx m — fit). 

The one-dimensional translational order of the smectic-A mesophase 
in addition to the orientational order found in nematics is characterized 
through the use of primary (orientational) and secondary (translational) 
order parameters together [50]. A complex order parameter can be used 
to characterize translational order [41]: 



where <j> is the phase, ip is the scalar amplitude of the density modu- 
lation, and A/B is the real/imaginary component of the complex order 
parameter. The density wave vector, which describes the average orien- 
tation of the smectic-A density modulation, is defined as a = V0/|V</>|. 
The smectic scalar order parameter ip characterizes the magnitude of the 
density modulation and is used in a dimensionless form in this work. In 
the smectic-A mesophase the preferred orientation of the wave vector is 
parallel to the average molecular orientational axis, n. 



2.2 MODEL FOR THE DIRECT ISOTROPIC /SMECTIC- 
A TRANSITION 



A two-order parameter Landau-de Gennes model for the first order isotropic/smectic- 
A phase transition is used that was initially presented by Mukherjee, 
Pleiner, and Brand [41,42] and later extended by adding nematic elastic 
terms [51,52]: 



f-fo= \a{Q:Q)-\b{Q-Q):Q + \c{Q:Qf 
+ |q|*| 2 + |/3|*| 4 
-±5|*| 2 (Q:Q)-feQ: (V*) (V**) 
+ \h (VQ) a + §ia (V-Q) 2 
+\hQ ■■ (VQ : VQ) 



where / is the free energy density, fo is the free energy density of the 
isotropic phase, terms 1-5 are the bulk contributions to the free energy, 




(1) 



* = ^ = A + iB 



(2) 



+ |6i|V*| 2 + ib 2 |V 2 4'| 
a = ao(T - Tjvi); a = a Q (T — Tai) 



(3) 
(4) 
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Figure 3: A schematic of one-dimensional direct isotropic/smectic-A growth 
under conditions where meta-stable nematic pre-ordering is present. The 
light/dark gray shading indication the magnitude of the nematic/smectic-A 
scalar order parameter S/ip growing in an unstable isotropic phase. Arrows 
indicate growth direction, where the meta-stable nematic front grows indepen- 
dently from the trailing smectic-A front. The increase in the nematic scalar 
order parameter is due to coupled enhancement from the presence of smectic 
order [42,53,55]. 



terms 6-7 are couplings of nematic and smectic order; both the bulk or- 
der and coupling of the nematic director and smectic density-wave vector, 
respectively. Terms 8-10/11-12 are the nematic/smectic elastic contri- 
butions to the free energy. T is temperature, Tni/Tai are the hypo- 
thetical second order transition temperatures for isotropic/nematic and 
isotropic/smectic-A mesophase transitions (refer to [53] for more detail), 
and the remaining constants are phenomenological parameters. This free 
energy density expression is a real-valued function of the complex order 
parameter ^ and its complex conjugate which makes it convenient 
to reformulate using the real and imaginary parts of the complex order 
parameter (see eqnj2|. 

As previously mentioned, the dual-order nature of the isotropic/smectic- 
A liquid crystal transition has been shown to result in meta-stable nematic 
pre-ordering. [3] shows a schematic of this phenomena, where a meta- 
stable nematic front precedes the stable smectic-A front, growing into an 
isotropic matrix phase. Under standard conditions, meta-stable nematic 
pre-ordering is not observed, but when the dynamic timescales between 
nematic/orientational and smectic- A/translational ordering differ greatly, 
meta-stable nematic pre-ordering is observed and predicted by the model 
based on eqn. [3] [37,39,43,54]. 

The previously mentioned dynamic timescales are included in this 
model via the Landau- Ginzburg time-dependent formulation [56]. The 
general form of the time-dependent formulation is as follows [56]: 




6 



where n„/fj, s is the rotational/smectic viscosity, and V the control volume. 
A higher order functional derivative must be used due to the second- 
derivative term in the free energy equation |3|: 

SF_ _ df_ _ _d_ ( _df_\ d d ( df \ 
56 ~ dd dxi I I + dxi dxj \ d 926 I 

where 6 corresponds to the order parameter. 



2.3 THERMAL ENERGY BALANCE DERIVA- 
TION 

The energy balance for a differential volume, without flow, is: 

^ = ^-V- 9 (8) 
dt dt w 

where u is the internal energy density, q is the heat flux, and is the rate 
of mechanical work on the differential control volume. According to eq. 
4.46 of ref [57] , the rate of mechanical work on the system can be written 
as the sum of elastic energy stored and the rate of dissipation. The rate of 
elastic energy storage is the change in nematic/smectic-A energy due to 
changes in the order parameters, so it is ^ taken at constant temperature. 
The rate of dissipation in the absence of flow is equal to [58]: 

dQ dA dB 

where A/B is the real/imaginary part of the complex order parameter 
(^), the molecular fields Hq, Ha, Hb are the negative of the variational 
derivatives of the free energy functional with respect to Q,A,B. With 
no flow, the Landau-de Gennes model gives (symmetric-traceless contri- 
butions are implied throughout this work): 

9Q df „ ( df 



dt dQ \8VQ 

The rate of mechanical work is thus: 

dW _ df dQ df ; dVQ 
dt dQ : dt + aVQ ' dt 

df dA df dT/A df 
+ dAl)t + dVA ' ~dT + : dt 

d£<m df dVB df dWB 
+ dB dt + dVB ' ~dT + 9VVB : dt 
, H . dQ dA dB 
+ HQ -^t +HA ^t +HB ^t 
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which simplifies to: 
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dt 



df dVQ 
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(12) 



Vavvs 

The left-hand side of eq[8]can be written in terms of Helmholtz free energy 

d 2 f 



du 
~dt 



d_ 

dt 



dTj 



d£ 
dt 



df_dT_ 
dT dt 



T 



dtdT 



(13) 



expanding the time derivatives yields: 



d£ 
dt 
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dAdT dt dBdT dt 



(14) 



(15) 



where it has been taken into account in eqn. [15] that all gradient terms 
in the free energy / are assumed temperature-independent. Additionally, 



the first term in eqn. 15 includes the specific heat (per unit volume): 



-T 



dT 2 



(16) 



Substituting eqns. [T4p6] into eqn. [13] yields the total change of internal 
energy: 



dt 



77 = Ci 



dT 

~dt 



T 



( d 2 f 
\dQdT 



dQ 

dt 



d 2 f dA d 2 f dB 



+ 



+ 



+ 



+ 



df dQ df ; dVQ 
dQ ' dt + dVQ' dt 
df_dA df dVA 
dA dt + dVA ' dt + 
df dB df dVB 



dAdT dt dBdT dt 



df 



dB dt + <9VB 



dt 



+ 



<9VVA ' dt 
df 9VVB 



<9VVB 



dt 



(17) 
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Finally, substituting the final expression for the total change in internal 
energy eqn. |17| and rate of mechanical work eqn. [12] into the energy 
balance eqn. [8] yields: 

9T 3Q 9Q (dA dA d_B dB_\ 

Cp dt ~ M ™ at : at + ^ s \dt : dt + dt : at ) 

T ( d 2 f dQ d 2 f dA d 2 f dB 
+ \8QdT ' dt OAdT dt OBdT dt 
- Vq (18) 

This can be further simplified taking the derivatives of free energy density 
eqn. [3] 

dT dQ 9Q (dA dA dB_ 8B\ 

Cp dt ~ Mn at : at + ^ s \dt : dt + dt : at ) 

/l _ dQ 1 ( .dA n dB\ 

- Vq (19) 

where term 1/2 corresponds to orientational/translational dissipation, 
term 3 corresponds to the energy of phase-ordering, and term 4 is the 
heat flux. A constitutive relationship for the heat flux that can be used 
is the anisotropic Fourier's law: 

q = -K ■ VT (20) 

where the thermal conductivity tensor K is used due to the anisotropy 
of the nematic phase. The thermal conductivity tensor can be written as 
the sum of isotropic and anisotropic contributions: 

/ fcii + 2k ± \ , , , , 

K = k lso S + k an Q = f 11 - ± J S + (fc|| - k±) Q (21) 

where ki SO and k arl are the isotropic and anisotropic contributions to the 
thermal conductivities, and fcii and k± are the conductivities in the direc- 
tions parallel and perpendicular to the average orientational axis (nematic 



director), respectively. While eqn 21 takes into account the dependence of 



thermal conductivity on the degree of orientational order [59,60], the effect 
of the translation ordering of the smectic-A phase on the thermal conduc- 
tivity is neglected based upon past experimental work. This work deter- 
mined that the shape anisotropy of the mesogens (liquid crystal molecules) 
is the main contribution to the thermal transport anisotropy compared 
to smectic layering [61]. Eqns. 3|5|18|21 constitute the model system of 
equations. 



2.4 SIMULATION CONDITIONS 

One-dimensional simulation of the model was performed using the Galerkin 
finite element method (Comsol Multiphysics). Quadratic Lagrange basis 
functions were used for the Q-tensor variables/temperature and quartic 
Hermite basis functions used for the complex order parameter compo- 
nents. Standard numerical techniques were utilized to ensure convergence 
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Figure 4: Schematic of the one-dimensional simulation domain where the initial 
smectic-A nucleus grows (from left to right) into the unstable isotropic phase. 
Both the density modulation (Re(^), solid line) and scalar nematic order pa- 
rameter (S, dotted line) are shown. The scalar smectic-A order parameter ip 
can be calculated using the both components of the complex order parameter 
t]j = VA 2 + B 2 



and stability of the solution including an adaptive backward-difference 
formula implicit time integration method. A uniform mesh was used such 
that there was a density of 6 nodes per equilibrium smectic-A layer in the 
liquid crystal computational domain of 500 layers (approximately 2/im). 
A coupled thermal domain, where only transient thermal diffusion was 
solved for, leveraged the fact that the characteristic thermal length scale 
(order of fim) is much greater than that of the characteristic smectic- 
A/meta-stable nematic length scale [62] (order of nm). Neumann bound- 
ary conditions were used at both boundaries and the initial condition 
consisted of uniform temperature with a single smectic-A layer nucleus. 
Additionally, a correction factor = 0.96 was used to normalize the slight 
difference between the latent heat value predicted by the phenomenolog- 



ical model parameters used (term 3 of eqn. 18 1 to an experimentally 
determined approximate reference value of AH exp — 4.85kJ/mol [63]. 
The latent heat predicted by the model and phenomenological parame- 
ters (determined previously [62]) is: 

df(T b ) „ /l _j, , 1 , 2 



AH = T„^^ = T b ^-a St + -aoiPt) (22) 

where Tb is the bulk isotropic/smectic-A transition temperature, based 
upon an approximate reference value experimentally determined [53] to 
be 331. 35^. 

Computational limitations constrain simulations to one-dimension in 
order to resolve the multiple length scales of smectic-A ordering (nm) and 
thermal diffusion (/im). Due to the use of the full complex order parameter 
(eq. [5| an additional constraint imposed on one-dimensional simulation 
requires that the smectic-A layers are parallel to the interface (imposed via 
initial/boundary conditions). A growing front where the smectic layers are 
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perpendicular to the interface is inherently two-dimensional and cannot 
be captured using one-dimensional simulation. 



3 RESULTS AND DISCUSSION 

Simulations were performed for three different cases: shallow isother- 
mal quench (neglecting the thermal energy balance, T — 331K), shal- 
low non-isothermal quench (T = 331K), and deep non-isothermal quench 
(T = 330K) where meta-stable nematic pre-ordering is predicted to be 
present (the timescales of nematic and smectic-A ordering differ greatly 
[43]). Shallow quench conditions describe the situation where the domain 
is quenched below the bulk isotropic/smectic-A transition temperature 
Tb but above the lower stability limit of the isotropic phase Tai- Sub- 
sequently, deep quench conditions describe a quench temperature below 
both T b and Tai- 

As previously mentioned, the vast majority of studies of diffusive dy- 
namics of liquid crystalline materials is focused on mass transport. How- 
ever, past experimental work by Dierking et al on growth laws of dif- 
ferent types of pure-component liquid crystal mesophases [33,34,64,65], 
including a study of the direct isotropic/smectic-A transition [34]. The 
culmination of this work, following the findings of Huisman and Fas- 
solino [32], includes a generalized non-isothermal model [66] showing the 
growth law evolution as a function of quench depth. This evolution 
transitions from purely non-conserved growth dynamics L oc t at deep 
quenches/undercooling to thermal diffusion-limited growth L oc t 1 ^ 2 un- 
der shallow quench/undercooling conditions. The presented extension eqn 
|18| of the high-order model of Mukherjee, Pleiner, and Brand [41,42] to 
account for latent heat effects allows for the resolution of the minimum 
physics |T]) to describe the direct isotropic/smectic-A transition ofa pure- 
component liquid crystal. 

[5^.-b shows the computed evolution of the power law exponent fit 
(I oc t n ) and interface velocities for the first two simulation cases (isother- 
mal and non-isothermal) versus position of the growing interface. [5^, shows 
that the inclusion of latent heat effects corrects the standard Landau-de 
Gennes model prediction of constant growth to that consistent with exper- 
imental observations of this specific system [34] and other liquid crystal 
mesophases [33,34,64,65]. [Sja shows that in addition to the correction 
to the growth law, the trend and magnitude of the interface velocity is 
modified. The initial relative increase in the front velocity is due to the 
fact that in this work nucleation effects are neglected and the initial nu- 
cleus size and degree of liquid crystalline ordering are assumed to be at 



bulk equilibrium values (see Section 2.4 1. As the stable smectic-A nucleus 



grows into the unstable isotropic phase, there is an initial decrease in the 
degree of smectic-A order as predicted by the high-order model. Due to 
the inclusion of latent heat effects this decrease in smectic-A ordering re- 
quires thermal energy, effectively resulting in local cooling. The interface 
velocity can be approximated by [2]: 

/3w = AF{T B -T*)-C (23) 
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Figure 5: Plots of the a) power law fit exponent (n in I oc t n ) and b) front 
velocity versus interface position for both the isothermal (solid line, T = 331 if 
and /is/^n = 25) and non-isothermal (dashed line, T = 331K and p s /fjt n = 25) 
simulations. The material parameters and phenomenological coefficients, based 
upon 12CB [62], are T NI = 322.85K, T AI = 330.5K, a Q = 2 x 10 5 



2.823 x 10 



7 J 



1.972 x 10 7 - 
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J,,a - 1.903 xl0 6 ^, 
5 = 9.792xl0 6 ^3, e = l.938xl6- n pN, l x = lxlCT 12 - 
bi = 1 x Kr 12 ^, 6 2 = 3.334 x 10- 30 Jm, and /x„ 
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2 pN,l 2 = 1.033xlO- 12 pAr, 
i>2 = o.ooi x ru - j in, anu /Lt n = 8.4 x 1CP 2 N ^ 2 S . Thermal 
material parameters used are based on values for 9CB and 10CB [67] (12CB data 
not available, T re f = 330. OK for thermal conductivity values) fcy = 0.3100- 
k ± = 0.1300-$*, C p = 2600^, and p = 1000^. 
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where /3 is the interfacial viscosity, AF (T s — T*) is the free energy driving 
force (free energy difference between the ordered/disordered phases) and 
C is the capillary force. T 3 and T* are the interface and transition tem- 
perature, respectively. Thus any local cooling/heating has a subsequent 
effect of increasing/decreasing the interface velocity. 

[6^,-b shows the evolution of the power law exponent fit (/ oc t n ) and in- 
terface velocities for the third simulation case (non-isothermal with meta- 
stable nematic pre-ordering) versus position of the growing interface. Un- 
der conditions where meta-stable nematic pre-ordering occurs, the effect 
of latent heat results in complex phase-ordering dynamics observed in 
other material systems [68]. 

In the present case, under deep quench/undercooling conditions, volume- 
driven growth is expected for the smectic-A phase. Conversely, any meta- 
stable nematic pre-ordering is under effectively shallow quench condi- 
tions (see Figure 2a of ref. [43] for the free energy landscape). As is 
seen in[6^, the smectic-A front initially undergoes constant volume-driven 
growth preceded by a meta-stable nematic front which rapidly exhibits 
diffusion-limited dynamics. As the meta-stable nematic front velocity de- 
creases |6]d) it is overtaken by the smectic-A front. This non-monotonic 
acceleration/deceleration/acceleration behavior is similar to dual phase- 
ordering/chemical-demixing growth kinetics recently simulated in nematic 
PDLC systems [68]. A unique differentiation in the growth dynamics 
as compared to those observed in the phase-ordering/chemical-demixing 
system is that the evolution of latent heat from the trailing smectic-A 
front contributes to the isotropic/metastable nematic interface tempera- 
ture which results in a small time period where the meta-stable nematic 
front begins to melt/recede (shown injfij}). 

[7] shows the interface temperature versus position of the growing inter- 
face for all three simulation cases. Inclusion of the latent heat contribu- 
tion from phase-ordering causes the interfacial temperature to approach 
the bulk transition temperature Tt, which subsequently causes the free 
energy driving force of the isotropic/smectic-A transition to approach 
zero. Thus, under these conditions, thermal diffusion limits growth re- 
sulting in dynamics converging to a power law exponent of n = | as 
observed experimentally [33,34,64,65] and predicted theoretically [32,35, 
66] for liquid crystal growth under shallow quench/undercooling condi- 
tions. In the case of metastable nematic pre-ordering, the evolution of the 
isotropic/metastable nematic front temperature reveals the influence of 
the latent heat contribution of the trailing smectic-A front. The thermal 
energy from the additional heat source causes non-monotonic behavior of 
the meta-nematic interface temperature evolution, favoring the eventual 
merging of the two fronts. Following this merging event, convergence of 
the isotropic/smectic-A front velocity is not accessible due to computa- 
tional limitations. The interface velocity is predicted to converge to a 
constant value due to the deep quench conditions and assumption of in- 
finite volume (see Section 2.4 1. These interesting dynamics suggest the 
possibility of oscillatory splitting/merging behavior of systems exhibiting 
metastable nematic preordering under realistic experimental conditions 
where ideal undercooling is not feasible. 
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Figure 6: Plots of the a) power law fit exponent (n in I oc t n ) and b) front 
velocity versus interface position for the non-isothermal simulation where meta- 
stable nematic pre-ordering is present (solid line for the smectic-A front and 
dashed line for the nematic front, T = 330K and n s /y>n = 250). The inset of 
b) shows the full scale of the nematic front velocity. 
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Figure 7: Plot of the interfacial temperature versus interface position of the 
three different simulations: non-isothermal (solid line), isothermal (dashed line), 
and non-isothermal with meta-stable nematic pre-ordering (stippled line for the 
smectic-A front and dash-dot line for the nematic front). 



4 CONCLUSION 

A thermal energy balance was derived to extend a high-order Landau-de 
Gennes type model of the isotropic/smectic-A liquid crystalline transition 
to take into account latent heat of phase-ordering, anisotropic thermal dif- 



fusion, and dissipation (eqn. 181. One-dimensional simulations were per- 
formed showing that the inclusion of the thermal energy balance corrects 
the standard Landau-de Gennes predictions of volume-driven growth un- 
der shallow quench conditions |5^i-b) to experimentally observed [33, 34] 
diffusion-limited growth. Additionally, the effect of incorporating ther- 
mal effects in the presence of meta-stable nematic pre-ordering was found 
to result in non-monotonic acceleration/deceleration/acceleration growth 
dynamics (Jg^l-Id) , as observed experimentally [37]. This work sets the basis 
for further study of non-isothermal effects of the direct isotropic/smectic- 
A transition where an extended subset of liquid crystal physics present 
in the transformation is included jl]). The presented extension to the 
high-order model of Mukherjee, Pleiner, and Brand [41, 42] allows for 
full three-dimensional multi-scale multi-transport simulation incorporat- 
ing the diverse defect and texturing dynamics of lamellar smectic-A or- 
dering coupled with non-negligible thermal and dissipative effects. 
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